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Abstract 

We consider a macroscopic body propagating in a one-dimensional resistive 
medium, consisting of an ideal gas at temperature T. For a whole family of 
collisions with varying degree of inelasticity, we find an exact expression for 
the effective force on the moving body as a function of the body's speed and 
the value of the restitution coefficient. At low and high speeds it reduces to 
the well-known Stoke's and Newton's law, respectively. 
Key words: Air drag, collisions 

Se considera un cuerpo macroscopico propagandose en un medio resistivo 
unidimensional, consistente de un gas ideal a temperature T. Para toda 
una familia de colisiones con diferente grado de inelasticidad, hallamos una 
expresion exacta para la fuerza efectiva sobre el cuerpo como funcion de 
la velocidad del cuerpo y del coeficiente de restitution. A bajas y altas 
velocidades, se reduce a la conocida ley de Stoke y Newton, respectivamente. 
Descriptores: roce viscoso, colisiones. 

PACS: 45.20.Dd , 45.50.Tn 
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1 Introduction 



When an object moves through a viscous medium, such as water or air, it 
experiences a resistive drag force. For small objects such as dust particles 
movingh at low speeds, this drag force is proportional to the speed of the 
object. This is known as Stoke's law|jl|. For larger objects such as airplanes, 
skidivers and baseballs moving at high speeds the drag force is approximately 
proportional to the square of the speed0. This limit is known as Newton's 
law. The general problem of determining the exact dependence of the drag 
force on the speed of an arbitrarily-shaped object moving at any speed, defies 
any closed form solution given its complex many-body character. A complete 
solution would have to take into account the detailed scattering process be- 
tween the body and the particles composing the medium, the thermal prop- 
erties of the medium, the presence of possible internal degrees of freedom 
of the body and local turbulence effects, etc. However, beneath all these 
complexities is basically the transfer of momentum and energy between the 
body and the medium particles. Therefore, is instructive to consider simpli- 
fied "toy" models where one can track in detail the momentum and energy 
exchange between the body and its surrounding medium. This is realized 
at the expense of simplifying other factors such as the dimensionality of the 
system or the specific form of the interactions between the body and the 
medium. In this spirit, we present here an extension of a previous 0, sim- 
plified zero-temperature one-dimensional model, where we now include finite 
temperature effects. This immediately brings into the game a natural veloc- 
ity scale not present in our previous model: the thermal speed. We obtain 
the resistive drag force as a function of the body's speed in closed form 
and find that, when the speed of the body is smaller than the thermal speed, 
the resistive force is linearly proportional to the speed of the body. On the 
contrary, when the speed of the body is greater then the thermal speed, the 
proportionality becomes quadratic. 
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2 The Model 



Let us consider a (macroscopic) body of mass M propagating in a one- 
dimensional resistive medium modelled by an ideal gas in thermodynamic 
equilibrium at temperature T, characterized by a thermal speed Vr = \J kT/m, 
where k is Boltzmann's constant and m is the mass of a medium particle 
(Fig. 1). We assume the body to be truly macroscopic, like a baseball mov- 
ing through air, or a falling rock. In other words, M ^> m which allows us 
to make the following simplification: During a medium particle-body colli- 
sion event, we will take the mass of the body to be essentially infinite. In 
this approximation the body is pictured as a massive, partially absorbing, 
moving "wall" colliding constantly with the medium particles. A reasonable 
assumption, if one considers that the mass ratio m/M is of the order of 10" 24 
for a baseball moving through air. After each collision, the speed of the body 
is essentially unchanged, so the magnitude of the momentum transferred to 
the body is 

Ap PS (1 + e)m\V — v\ 

where V is the speed of the body, v the speed of the medium particle and e 
is the restitution coefficient for the body-particle collision. Thus, when e = 1 
we have a completely elastic collision, where the magnitude of the relative 
body-particle velocity is conserved, while at e = 0, we have the case of a 
completely inelastic collision, where the particle is "absorbed" by the body 
after colliding. We also work in a quasi-continuum approximation where an 
element of length dx while "small" , will contain a large number of medium 
particles. 

Initially the body is given a speed V (say, to the right), and we observe 
the system at a later time t, when the speed of the body is V. During the 
next time interval dt, the body will collide with particles coming from its 
left and right side. On the left side, only those particles that have speeds 
v > V and are located closer than (v — V)dt will collide with the body. The 
number of such particles is driL — pdn(v)Q(v — V)(v — V)dt. In a similar 
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manner, the number of particles to the right of the body that will collide 
with the body during the interval dt is driR = pdn(v)Q(V — v)(V — v)dt. 
Here p is the particle number density, Q(x) is the step function (O(x) = 1, 
for x > 0, zero otherwise) and dn(v) is the number of particles that have 
speeds in the interval [v,v + dv\. dn(v) = g(v)dv where g(v) is the thermal 
speed distribution, given by 

1 



[l/V2n)-eM-\(v/V T ) 2 } 



//(./') ■ (i/ V^ttj— cxp[--(/v» r ) |. (1) 

The transfer of momentum per unit time coming from the medium to the 
left of the body, due to particles with speed in the interval [v, v + dv] would 
then be: 

dP 

— = (1 + e)m(v - V)g{v)e{v - V)(v - V)dv (2) 

By integrating over all speeds, we obtain the average effective force on the 
body from the left side: 



left 



pm(l + e)(v — V) 2 g{v)dv. 



(3) 



In a similar manner, the transfer of momentum per unit time coming from the 
medium to the right of the body, due to particles with speed in the interval 
[v , v + dv] is: 



dP 

_ = -(! + e )m(V - v)g(v)Q(V - v)(V - v)dv 



(4) 



which implies that the average force on the body from the right side is 

bright = [ p(l + e)(V-v) 2 g(v)dv (5) 



The net average force F on the body, along the direction of its initial velocity, 
is given by the difference between Eq.(|3|) and (|5|): 

rV 



-mp(l+e) 



/V poo 
p(l + e)(V - vfg{v)dv - / pm(l + e)(v - V) 2 g{v)dv 
-co JV 



By inserting expression ([l|) for g(v) and carrying out the integrations, we 
obtain: 



-mp(l+e)V q 




1 /V- 1 
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1+i B) Eli {V2V T 
(6) 



This rather complex-looking expression is a bit deceiving since it depends on 
negative exponentials of (V/Vr) 2 which makes it very sensitive to whether 
V/Vt < 1 or V/Vt > 1. In other words, we expect two, well-defined behavior 
regimes, with a small crossover region near V/Vt ~ !■ 

3 Results and discussion 

As can be clearly seen from (P), the degree of inelasticity plays only a minor 
role, renormalizing the number density of the medium. Figure 2 is a log-log 
plot of the effective average force on the body as a function of the speed of the 
body, Eq.(||). As anticipated above, we note that except for a small vicinity 
around V = Vr, it consists of basically two straight lines with slopes of one 
and two, respectively. That is, at speeds smaller than the thermal speed Vr, 
the resistive force is proportional to the body's speed (Stoke's law); while for 
body's speeds greater than Vr, the resistive force becomes quadratic on the 
body's speed (Newton's law). These limits are easy to derive from Eq.(|6]): 
For V < V T , Eri(V/V2V T ) « y/2/n(V/y/2V T ) and exp[-(l/2)(V/V T ) 2 } w 1, 
which implies: 

F » -yj8/Trmp(l + e)V T V = -p(l + e)^8mkT/n V V < V T . (7) 

On the other hand, when V > V T , Erf(V/y/2V T ) ~ 1 and exp[-(l/2)(V/V T ) 2 ] « 
0. Thus, in this case one has: 

F « -mp(l + e)V 2 V > V T . (8) 

Let us now consider the issue of the stopping distance. For a medium at a 
finite temperature, the speed of the body decreases (on a macroscopic scale) 
as it moves through the medium and will eventually become smaller than 
the thermal speed. At that point, the resistive force becomes proportional 
to the speed, F = —f3 V. A simple integration then leads to an exponential 
decrease on V and therefore, a finite stopping distance. If the medium is at 
zero temperature however, the resistive force is always quadratic with speed 
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F = jV 2 and, in that case, it can be easily proved that the stopping distance 
diverges logarithmically with time||. 

In summary, we have examined a simplified model of a macroscopic object 
propagating in a resistive one-dimensional medium modelled as an ideal gas 
at temperature T. For general inelastic collisions between the body and the 
medium particles, characterized by a restitution coefficient e, < e < 1, we 
have arrived at a closed-form solution for the resistive force in terms of the 
speed of the body. Below the thermal speed, this force is essentially linear 
in the body'd speed, while above thermal speed, the dependence becomes 
quadratic. 
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Figure Captions 



FIG 1: A macroscopic body of mass M propagating inside a one-dimensional 
resistive medium composed of an ideal gas of particles of mass m, with m <C 
M, in thermal equilibrium at temperature T. The body undergoes partially 
elastic collisions with the medium particles with restitution coefficient e. 

Fig. 2: Effective average force on the macroscopic body as a function of 
the body's speed. For speeds smaller (higher) than the thermal speed, the 
dependence is essentially linear (quadratic). The crossover region is confined 
to a small vicinity around Vt- ( F = mp{l + e)V^ ). 
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http://arXiv.org/ps/physics/0207037vl 




FIG. 2 



